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In this review, three subjects are discussed: a phenomenological application of lattice predictions to K^-K'^ mixing 
in Super Symmetry; a discussion of the non-perturbative renormalization methods for four-fermion operators and 
a new approach to extract weak matrix elements in effective theories denoted as OPE without OPE (operator 
product expansion without operator product expansion). 
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1. Introduction 

Since the original proposals of using lattice QCD 
to study hadronic weak decays substantial 
theoretical and experimental progress has been 
made: with both Wilson and staggered Fermions, 
the main theoretical aspects of the renormaliza- 
tion of composite four-fermion operators are, by 
now, fully understood the calculation of the 
K^-K'^ mixing amplitude, expressed in terms of 
the so-called renormalization group invariant B- 
parameter Bk^ has reached a level of precision 
(at least in the quenched approximation) which 
is unpaired by any other approach @-||] (for 
a detailed discussion see [^); increasing preci- 
sion has also been gained in the determination 
of the electro-weak penguin amplitudes (and of 
the strange quark mass JlOt), which are relevant 
in the prediction of the CP violating parameter 
e'/e. A lot of progress has also been made in 
the determination of quantities which enter heavy 
hadron weak processes and which are reviewed by 
T. Draper in his plenary talk Still, in par- 

ticular for kaon physics, many problems are un- 
solved or poorly understood. By using chiral per- 
turbation theory and SU{2i) symmetry, for exam- 
ple, one may relate Bk to the physical A/ = 3/2 
K — > TTTT amplitude, A2. The large value of 
Bk ^ 0.85 found in lattice calculations leads to a 
prediction for A2 which is about 40% larger than 
its experimental value. It is not clear whether this 
is due to a failure of chiral perturbation theory, 
or to electro- magnetic effects |]l2|. In the former 
case, the use of X ^ vr matrix elements, com- 
bined with chiral perturbation theory, to predict 



the A/ =1/2 amplitude or e'/e would be very 
suspicious. Until very recently, the only exist- 
ing calculation of two other i?-parameters which 
are very important for e'/e, namely B^ and Bg, 
was the one performed with staggered Fermions 
in iQ. Then, the perturbative one-loop lattice 
corrections to the relevant operators were com- 
puted and found so large as to make the old 
estimate of Q very doubtful (a recent calcula- 
tion of was presented in |15|). Thus, we do 



not have any reliable lattice estimate of B^ and 
Bq. Finally, there is no convincing lattice calcu- 
lation of Aq yet, although some new results with 
staggered fermions have been presented at this 
Conference @. 

In spite of considerable technical improvements, 
and of the increasing computer power, no fun- 
damental advance has been done last years for 
the problems mentioned above. For this reason, 
rather than presenting new results, I discuss new 
methods which have been proposed in order to 
compute the relevant matrix elements. In most 
of the cases, these methods have not been applied 
yet, or only feasibility studies exist. This is the 
content of sec. ^. In sec. ||, I present a recent ap- 
plication of lattice calculations of AS* — 2 ampli- 
tudes to flavor changing neutral current (FCNC) 
phenomenology in Super Symmetric models. In 
sec. ^ several approaches, which have been used 
or proposed for the non-perturbative renormal- 
ization of lattice four-fermion operators, are dis- 
cussed and compared. 
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2. AS = 2 transitions in SUSY 

The prescription of minimality in tlie number 
of new particles introduced to supersymmetrise 
the Standard Model (SM), together with the de- 
mand of conservation of baryon and lepton num- 
bers, does not prevent the appearance of more 
than 100 new SUSY parameters, in addition to 
the 18 aheady present in the SM. FCNC and 
CP violating phenomena are protagonists of a 
drastic reduction of these extra degrees of free- 
dom. Among the possible FCNC processes, K'^- 
j^o (^bO_j^o^ mixing play a very special role. In 
this section, a recent analysis of K^-K^ mixing in 
SUSY is discussed jl^. This analysis makes use 
of the recently computed NLO corrections to the 
SUSY effective Hamiltonian and of the ma- 
trix elements of the relevant operators (renormal- 
ized non-perturbatively in the RI-MOM scheme), 
computed on the lattice . This is a new, inter- 
esting application of the lattice method for study- 
ing the physics beyond the SM. Unfortunately, for 
lack of space, many other applications to physics 
beyond the SM cannot be discussed here. A de- 
tailed discussion of the non-perturbative meth- 
ods which can be used to renormalize the relevant 
four-fermion operators can be found in sec. ^. 
Let us start by defining the model used to con- 
struct the low-energy effective Hamiltonian for 
AF = 2 transitions. In the so-called mass inser- 
tion approximation |^^, one chooses the super- 
CKM basis for the fermion and sfermion states, 
where all the couplings of these particles to neu- 
tral gauginos are flavour diagonal, while the gen- 
uine SUSY FC effects are exhibited by the non- 
diagonality of the sfermion mass matrices. Denot- 
ing by the off-diagonal terms in the sfermion 
mass matrices (i.e. the mass terms relating 
sfermions of the same electric charge, but differ- 
ent flavour), the sfermion propagators can be ex- 
panded as a series in terms oi 5 = A? /rh? , where 
TO is the average sfermion mass. As long as A^ 
is significantly smaller than fh? , we can just take 
the first term of this expansion and, then, the 
experimental information concerning FCNC and 
CP violating phenomena translates into upper 
bounds on the 5s [^~|22|. There exist four differ- 
ent A mass-insertions connecting flavours d and s 
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Figure 1. Feynman diagrams for A5 = 2 transi- 
tions, with h, k,l,m = {L, R}. 



along a sfermion propagator: (Afj)^^, (Af2)^^, 
(Afa)^^ and (Afa)^^- The indices L and it! re- 
fer to the helicity of the fermion partners. For 
gluino-mediated processes, for example, the am- 
plitude for AS* = 2 transitions in the full theory 
is given by the computation of the diagrams in 
flg. |l|. Having calculated the amplitude from the 
diagrams in fig. |l|, one has to choose a basis of lo- 
cal operators and perform the matching of the full 
theory to the one described by Hf/=2. In [0|9|, 
they have adopted the form 



AS=2 
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and the operators Qi,2,3 are obtained from the 
Qi,2,3 by the exchange L ^ R. Here qj^ l = 
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PR,Lq, with Pr^l = (1 ± 75)/2, and a and (3 
are colour indices. The Wilson coefficients Ci(/i) 
(Ci{iJ,)) for the operators renormalized at the 
scale n are obtained by computing the coefficients 
at the large energy scale, Ci{Ms), where Ms is 
the order of the gluino (squark) mass, and then 
evolving them to /Lt using suitable QCD renor- 
malization group equations, which depend on the 
operator anomalous-dimension matrix. The lat- 
ter is only determined from QCD and does not 
know about the underlying fundamental theory 
at the large energy scale. The next-to-leading 
anomalous- dimension matrix for the most gen- 
eral Ti^^^ has been computed in Q. They 
used the regularisation-independent (RI) anoma- 
lous dimension in the Landau gauge, since we will 
make use of matrix elements computed in lat- 
tice QCD with the same choice of renormalisa- 



tion scheme 19 1. A full NLO computation would 
also require the 0{as) corrections to the match- 
ing conditions. Unfortunately, such corrections 
are not available yet. One might argue that, be- 
ing of order as{Ms), these contributions should 
be small, as suggested by the cases of the SM and 
of the two Higgs doublet model; this statement 
can only be confirmed by an explicit computa- 
tion. Due to the absence of 0{as) corrections 
to the matching, our H^^"^ will be affected by 
a residual scheme dependence, which would be 
cancelled by the missing terms of order as{Ms). 
A simple analytic formula for the expression of 
the Wilson coefficients function of the ini- 
tial conditions at the SUSY scale, C(Ms), and of 
as{Ms) can be given. For Ms > rrit, one has 

= E E (^^"^ + 

i s 

where, in the evolution of the coefficients from 
Ms, Ms = {Mg + Mq)/2 has been chosen. 77 = 
as{Ms) / as(jnt), and the magic numbers a^, b\'^ 
and c^'* at /i = 2 GeV in the RI-MOM scheme can 
be found in [pT| . The magic numbers for the evo- 
lution of Ci_3 are the same as the ones for the 
evolution of 6*1-3. Eq- (§) can be used in con- 
nection with the i?-parameters, given in eq. (|^) 
below, to compute the contribution to /S.Mk and 
Ek at the NLO for any model of new physics in 
which the new contributions with respect to the 



SM originate from extra heavy particles. One has 
just to plug in (H) the Ci of his favourite model. 
In the case of the renormalised operators, a con- 
venient definition of the _B-parameters is the fol- 
lowing 
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where the notation Qi{fi) (or simply Qi) denotes 
the operators renormalised at the scale /i. Here 
and in the following, the same expressions of the 
i?-parameters of the operators Q1-3 are valid for 
the operators Qis- In eq. (||) the operators and 
the quark masses are renormalised in the same 
scheme (RI, MS, etc.) at the scale /i and the nu- 
merical results for the i?-parameters, Bi — Bi(fj,), 
presented below refer to the Landau RI scheme. 
Moreover, without loss of generality, terms which 
are of higher order in the chiral expansion and 
which are usually included in the definition of the 
i?-parameters have been omitted. The advantage 
with the definition (jj) is that the S-parameters 
obeys to very simple renormalization group equa- 
tions, with an anomalous dimension matrix which 
is related in a simple way to that of the corre- 
sponding operators |p^ . In the numerical study 
below, the following values of the i?-parameters, 
for /i = 2 GeV, have been used 



Bi = 0.60(6) B2 = 0.66(4) 

B3 = 1.05(12) 

Bi = 1.03(6) B5 = 0.73(10) 



(5) 



The central value used for Bi — Bk corresponds 
to Bf(^^{2GeV) = 0.61, in agreement with the 
recent estimates of B2-5 have been taken 
from lis] , where all details of the computation can 
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be found (for another determination of these B- 
parameters, calculated with a perturbative renor- 
malization, see [p3|). 

As an example, I present the results of a model- 
independent analysis of i^^-J^^ mixing. The Kl- 
Ks mass difference AMk and the CP-violating 
parameter sk are given by 

AMk = 2Re {K'>\n^s^=^\K°) , (6) 

The SUSY (gluino-mediated) contribution to the 
low-energy H^ff^^ contains two real and four 
complex unknown parameters: rriq, rrig, {^12} j^j^, 

{6(2) LW (OflL and {Si^)^^. The model- 
independent constraints are obtained by imposing 
that the sum of the SUSY contributions propor- 
tional to a single (5-parameter and of the SM con- 
tributions to AMk and ek does not exceed the 
experimental value for these quantities. This is 
justified by noting that the constraints on differ- 
ent ^-parameters in the kaon case exhibit a hierar- 
chical structure, and therefore interference effects 
between different contributions would require a 
large amount of fine tuning. The B-parameters 
used in the analysis are those given in eq. (j^) sub- 
tracted by one standard deviation to their cen- 
tral values, in order to extract a more conserva- 
tive bound on SUSY parameters. The limits on 
the 5-parameters for a typical value of the gluino 
and squark masses are reported in tables and ^ 
(more detail can be found in ||l7[]) . 
As a glimpse at tables |l| andg readily shows, 
the inclusion of the NLO QCD corrections to the 
Wilson coefficients of H^g^^ and the use of the 
lattice i3-parameters, instead than the matrix el- 
ements evaluated in the vacuum saturation ap- 
proximation (VSA) , affects the results in different 
ways, according to the different operators. The 
effects are particularly large for left-right opera- 
tors. 

3. Non-perturbative renormalization of 
four-fermion operators 

From the theoretical point of view, the calcula- 
tion of hadronic weak decays can usually be di- 
vided in two parts: the calculation of the effects of 
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Table 1 

Limits on Re((5y)^^ i^ij)cD' '^ith A,B,C,D = 
{L,R), for an average squark mass rriq = 500 GeV 
and for different values of a; = m^/nr-. 
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Table 2 

Limits on Im((5y)^^ i^ij)cD^ with A,B,C,D = 
(L, i?), for an average squark mass niq = 500 GeV 
and for different values of x = rn^/ra^. 
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strong interactions at short distances, which can 
be done in perturbation theory using the Wilson 
Operator Product Expansion (OPE), and the cal- 
culation of the hadronic matrix elements of local 
operators. The latter contain the effects of strong 
forces in the infrared region where perturbation 
theory cannot be applied and the lattice approach 
is used. This is the path followed in sec. I for the 
costruction of H^ff^^. An important point in 
the full approach is the consistency between the 
procedure used to obtain the renormalized opera- 
tors and the scheme adopted in the calculation of 
the Wilson coefficients. Unlike in many other for- 
mulations, as for example in the 1/N expansion 
beyond the leading order [Q, and in the chiral 
quark model p5[ , the matching of the renormal- 
ized operators to the corresponding Wilson coef- 
ficients can be rigorously done, at least in princi- 
ple, using lattice perturbation theory. In practice, 
the matrix elements of the renormalized operators 
constructed from the bare lattice ones are sub- 
ject to two main sources of systematic effects: the 
renormalization constants (matrices) are usually 
computed in one-loop lattice perturbation theory 
and are subject to large C(g^) errors; in most of 
the cases the matrix elements are computed us- 
ing the Wilson action and they suffer from 0{a) 
discretization errors. 

Numerical studies of the matrix element of 
the AS = 2 left-left operator, 0^-5=2 ^ ^^^^^ - 
75)^57^(1 — 75)c?, have shown that the chiral be- 
haviour of the renormalized operator is not signif- 
icantly better if one adopt a tree-level improved 
action and operators, whereas it is significantly 
improved by using non-perturbatively determined 
renormalization constants (either by Ward iden- 
tities [ p6[ or with the non-perturbative method 
(NPM) of |2|)[]. 

In this section, I review some of the methods 
which is possible to use in order to renormal- 
ize non-perturbatively the lattice operators, thus 
eliminating one of the main sources of uncertain- 
ties in the calculation of the physical amplitudes. 
To this purpose it is useful to start with some 
definition and classification of the relevant four 



^ For a different conclusion, using a boosted improved ac- 
tion and operators, see L. Lellouch at this Conference. 



fermion operators. In the following I only discuss 
the case of AS* = 2 or AI = 3/2 operators, which 
cannot mix with operators of lower dimensions. A 
full discussion of the renormalization of A J =1/2 
operators can be found in and more recently 
in 0, see also 

Operator basis In order to illustrate the renor- 
malization properties of these operators, it is con- 
venient to work in a theory with four distinct 
flavours and define the generic operator as 



O 



rir2 = \ [OtiF: 



(8) 



2 [(■0iri'02)(V'3r2'04) ± {lpl^l'>p4){lp3^2tp2)] 

where ri^2 denote one of the 16 Dirac matrices, 
r = {I,7M''^M'''7/i75,75} = {S,V,T,A,P}, and 
summation of Lorentz indices is implied when 
necessary. For example Ovv corresponds to 
yV = E^7/x 7/x, VA = E^7m 7m75, etc. 
In total we can classify ten parity-even and ten 
parity-odd operators. There are several possible 
choices for the operator basis. In the following, I 
will use 



Qf 



± _ 



'-^[VV+AA] 
^[SS-PP] 



Q2 — ^[VV~AA] 
^4 — '^[SS+PP] 



Qf = a 



TT 



(9) 



for the parity-even operators and 



^[VA+AV] 
^[SP-PS] 



^2 — '^[VA-AV] 
^^4 — ^[SP+PS] 



(10) 



for the parity-odd ones. 

Parity, chirality and generalized CPS symmetries 
restrict the possible mixing between different op- 
erators [|l|. The simplest case is that of the 
parity-odd operators for which the mixing ma- 
trix, Z.y , is a sparse block diagonal one. We have 



(11) 



where the Q,; denote the renormalized operators 
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and 
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(12) 



The block-structure of the matrix Zij is not 
changed by a lattice regularization which breaks 
chiral symmetry. The situation is rather differ- 
ent for the parity-even operators Qi. In an ideal, 
chirally symmetric regularization scheme, xi?S', 
the mixing matrix of the parity-even operators 
Zij would have the same block-structure as the 
mixing matrix of the parity-odd ones. The ex- 
plicit chiral symmetry breaking induces further 
mixing which can be removed by suitable counter- 
terms which enforce, up to discretization errors, 
the relevant Ward identities for the subtracted 
operators. A convenient parametrization of the 
renormalization matrix is given by the following 
expressions 



= [/ + A±]Q± 



(13) 
(14) 



where 
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and 
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) 





(15) 



(16) 



Note that in the hypothetical xRS, A^ = and 

Ward identities and the NPM I now dis- 
cuss, with one simple example, the determination 
of the mixing matrix A^ using the Ward iden- 
tity method (WIM). The WIM, first discussed 
in Q (see also Q), was apphed in Q to de- 
termine the renormalization constant of the axial 



current using quark Green functions and then im- 
plemented to the non-perturbative determination 
of the mixing coefficients of the AS* = 2 opera- 
tor Qi in (2|]. I wiU also show that the WIM is 
equivalent to the NPM on quark and gluon states, 
first proposed in [2^ and then implemented for 
the renormalization of the four-fermion operators 
in |27|. A more detailed discussion can be found 
in 1 34]. Before discussing the Ward identities, I 
have to introduce some definition. The Green 
functions considered in the following are given by 
the expectation values of the following multilocal 
operators: 



Gfc(xo; 2:1,2,3,4) 

Qk{xo]Xi^2.,?.A) 

tjf ^(xo;a;i,2,3,4) 
^^''(xo; 2:1,2,3,4) 



'0iV'2Qfe(a;o)'03V'4 

i'i'4^2 Qkixo)i>3i>4: 

i>l[lp2l5]Qkixo)lp3lp4 



and similarly for Q^^"^ and Gf^^ k — 1, ... ,5. In 
the above equation tpi = ^pl{xl), 1^2 = "02(2:2), 
etc. For simplicity, the colour and spin indices 
have not been shown explicitly and the ± su- 
perscripts have been dropped from operators and 
correlation functions. Having thus defined Gfc, Qk 
etc. in coordinate space, we will also be using 
the corresponding Green functions in momentum 

space Gfe(p), Gk{p), Gk^\p) etc. (all external legs 
having the same Euclidean momentum) . The am- 
putated Green functions of {Gk{p)), {Gk{p)) and 
{G'k\p)) ^ill be denoted by Afe(p), 3\i.{p) and 
3\^k\p)^ respectively. 

On quark states, the relevant Ward identity for 
the operator Qi can be written as 



d^x{Gi [V^A^(:e) - 2moP(x) - X^{x)]) 



(17) 



where the flavour indices of Afj^{x) and P{x) 
have been omitted. Out of the chiral limit, the 
2moP(a;) term is present, but the surface term 
from the current divergence vanishes upon inte- 
gration. I will be considering this case, in or- 
der to mimic what is practicaly done in numeri- 
cal simulations (i.e. first we compute quantities 



7 



at small non-zero quark mass and then extrapo- 
late to the chiral limit). The operator X5 arises 
from the variation of the chiral symmetry break- 
ing Wilson term in the action. As shown in 
(see also for a detailed discussion) it mixes, 
under renormalization with V^Ap, and P. This 
mixing, determined by the requirement that on- 
shell matrix elements of the subtracted X5 vanish 
in the continuum, generates a finite renormaliza- 
tion of the axial current and a power subtraction 
of the quark mass. Thus, following in the 
above WI we will trade-off X5 for the renormal- 
ized expression [V^A^ — 2mP — X5], where 
is the subtracted X^. What is of interest to us 
is that, besides the above renormalizations, the 
X5 insertion in the above corelation function also 
generates new contact terms 



(l-|ii)(Gi)<5(:r-xo) 
^11 

-|ii^Aifc(Gfe)J(a;-xo) 



(18) 



where the notation for the various coefficients 
has been chosen with some foresight. There are 
also contact terms arising from the proximity of 
Xc^(x) to the quark fields of the correlation Qi 
at points xi, ■ ■ ■ ,Xi. These terms have the form 
{Q^'')5{x — Xk) (with k — 1, • • • ,4). However, as 
shown in | ]36[ |, they vanish in the continuum limit. 
We now combine eqs. (0) and (|l8t), Fourier- 
transform the WI (with all external momenta set 
equal to p) and amputate the resulting correlation 
functions 



Ai(p) = Ai(p) + AifeAife(p) 



k=2 



4Zn 
4' 



'-^ d^x{g^{p)mP{x))\{{S-\p)) (19) 
11 -> j=i 



We require the above WI, up to quark field renor- 
malization, to be the one valid in the continuum 
limit for subtracted correlation functions (oper- 
ators). This implies that the l.h.s., multiplied 



by Zii, can be identified with the renormalized 
parity-conserving correlation function (operator). 
We see immediately that the WI fixes the opera- 
tor subtractions (i.e. the mixing coefficients Ai^). 
This can be done by projecting the above WI 
with suitable projectors |Q and by solving the 
resulting system of linear non-homogenous equa- 
tions for the four A's ^ The WI also fixes the 
ratio of the multiplicative renormalization con- 
stants Zii/Zii, but not each of them separately. 
This is the analog of what happens in the case of 
the scalar and pseudoscalar densities: the Ward 
identities can only fix the ratio of the renormal- 
ization constants of these operators, Zp/Zs, but 
not Zp or Zs- Their values, for fixed Zp/Zs, 
is arbitrary and depend on the renormalization 
scheme. 

The WI (|l^) is only a function of bare quanti- 
ties. This demonstrates that the matrix A and 
the ratio ZwjZw are finite functions of ^q, which 
can be completely determined from the Wis, and 
do not depend on the renormalization conditions. 
The overall renormalization constant Z\\ (-Zn), 
instead, is logarithmically divergent (as in the 
continuum) and is fixed by imposing a suitable 
(scheme-dependent) renormalization condition on 
the subtracted operator Q\. 

I will demonstrate now that the WIM and the 
NPM are equivalent for the determination of the 
lattice mixing coefficients. The validity of this 
statement has already been discussed in general 
in and in the specific case of the four-fermion 
operators in ||3^. Here I discuss the specific ex- 
ample of Q\. It is to be understood, throughout 
the rest of this section, that the chiral limit is to 
be taken in the end. 

From (p^, in the large momentum limit, one may 
derive the renormalization conditions used in the 
NPM to fix the matrix A. For large p, the last 
term on the r.h.s. of eq. (|l9|) is power supressed. 
This happens because the explicit m factor im- 
plies that the integrand has one less dimension 
that the other terms so that, at large momenta. 



The projectors obey the following orthogonality condi- 



tions: Tr Pf A 



±A±(0) 



5,fc and Tr P±>AJ'°' 



j,A; = 1, 



. , 5 where A^°>=^ and 



with 
are the tree-level 



amputated Green functions of q¥ and respectively. 
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it vanishes faster by one power of p. Moreover, 
in this hmit the inverse quark propagator has the 
general form S~^{p) = iSi7^p^ (with Ei a scalar 
form factor). This means that in this limit it an- 
ticommutes with 75. By projecting both sides 
of eq. ( p^ ) on the single operators Qi, (with 
k = 1, . . . ,5) one obtains 



Tr P, Ai 



5 

E 

k=l 



AuTr PjAk 



= |^Tr 1-^1(20) 
^11 



which is precisely the system of equations used 
in and to fix the coefficients of the mixing 
matrix A. 

The generalization to the other operators is 
straighforward, although more complicated. For 
example, considering the operators Q2.3 and Q2,3, 
we have 



[A(p) + Aa(p)zI]zJ = --[:A^'^ 



4 

d^x{g{p)C^2mP{x))l[{S-\p)) (21) 



where I have used a compact notation for the 
row- vectors G = (^2,^3), A = (A2,A3), = 
iJ\2,J\3), Aa = (Ai,A4,A5) and for the mixing 
matrices 



ZA 



A21 
A31 



^22 


^23 


Z32 


Z33 


Z22 


Z23 


Z32 


Z33 


A24 


A25 


A34 


A35 



(22) 



Equation (|21j) correspond to a system of linear 
equations in the matrix elements of za and in 
^ which can be used to fix the subtracted 
operators Q2,3- Finally, overall renormalization 
conditions have to be imposed to the Q2,3 to ob- 
tain the mixing matrix z^. 
Finally, I want to mention another method, im- 
plicitly suggested in which has not been 
implemented yet in numerical simulations. This 
method is based on the Ward identities which 
can be written for gauge-invariant correlation 
functions. Let us introduce the "meson fields" 
P^itx) - j (fxs{x,tk)i5d{x,tk), P^'+iU) = 



J (Pxu{x,tTr)^5d{x,tjr), ctc. By defining the fol- 
lowing correlation functions 

G2{U„U,) = {0\P'''{U,)Q{xo)P''HU,)\0) 

G3{U,,U,)^ J dtK{0\P''{tK)Q{xn) 

P'^'iUjP'^'iU.m (23) 

the Ward identity can schematically be written 



lim 2TOZtJ3(tTri, ^TTa) = const. X ZG2 (^Tri, ^TTa) (24) 

m — >0 

where const, is a suitable constant which depends 
on the flavor quantum numbers of the operators Q 
and Q and of the external "meson fields" . For dif- 
ferent values of Itti and , the above Ward iden- 
tity results in a system of linear non-homogenous 
equations from which it is possible to extract A='= 
and ZyZ^ ^ . A study in this direction is under- 
way 

4. OPE v^fithout OPE 

In this section, I will discuss some new methods 
which have been developed in order to overcome 
the difficulties related to the renormalization of 
lattice composite operators relevant in many im- 
portant physical applications. Among the others, 
let me mention the structure functions in deep in- 
elastic scattering, kaon and i?-meson non-leptonic 
exclusive decays, electromagnetic form factors at 
large momentum transfer, exclusive and inclu- 
sive semileptonic i?-decays etc. The general idea 
underlying these new methods is essentially the 
same and I will denote it as Operator Product 
Expansion without Operator Product Expansion 
(OPE without OPE or OPEwOPE). The reason 
for choosing this name will be clearer in the fol- 
lowing. I now discuss the general features of the 
OPEwOPE, starting from the description of the 
standard approach followed so far in lattice cal- 
culations. 

Let us consider the short-distance OPE for the 
T-product of two generic currents between the 
physical external states \A) and \B) 

{B\T{Jix)jmA) ^ ^a(a:,/.)(B|0.(M)|A) (25) 
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where, for simplicity, Lorentz indices have been 
omitted. In the above equation, ^ is the renor- 
mahzation scale of the local, composite operators 
Oi(/Lt) and the Ci{x, /i) are the Wilson coefficients, 
which can be computed in perturbation theory, 
provided fi ^qcd- The Wilson coefficients 
are choosen in such a way that the T-product is 
fi independent. The standard approach consists 
in computing the Wilson coefficients in perturba- 
tion theory and extracting the matrix elements of 
the local operators from suitable hadronic corre- 
lation functions. Many examples of matrix ele- 
ments of local operators that have been studied 
in the past can be given: i) Moments of the struc- 
ture functions, e.g. {pl'ipij^ii Df_,2 ■ ■ - i i/'b); 
ii) Hadronic weak decays {K tttt or _B — > tttt), 
e.g. (ttttI (s7^d) (uj^e) \K); in) K°-K° and B°- 
B° mixing, e.g. {R'^l (sjf^d) (sj^d) {K'^); iv) 
operators of the HQET entering inclusive and 
exclusive decay rates, e.g. {B\hyDohy\B) and 

{B\hy{i D)^hy\B); iv) light-cone wave functions 
for electromagnetic form factors and exclusive 
B TT and B ^ p semileptonic decays, e.g. 

The renormalized operators Oi{^) are obtained 
from the bare lattice operators either by us- 
ing (boosted-tadpole improved) perturbation the- 
ory or non-perturbatively using the WIM or the 
NPM. The difficulties encountered in this ap- 
proach are due to the fact that the ultraviolet be- 
haviour of the "effective theory" (OPE, HQET, 
NRQCD, etc.) is much worse and divergences 
which were not present in the "full" theory ap- 
pear in the different expansions. For example, if 
the currents are the SM quark currents responsi- 
ble for weak decays, the l.h.s. of eq. ( p5|) is finite, 
whereas the operators appearing on the r.h.s. are 
logarithmicaly or power divergent in the inverse 
lattice spacing a^^ . We briefly discuss examples 
taken from the cases considered above, 
i) Given the bare energy-momentum tensor den- 
sity T^'' = '0(7^ Du +lu D^)ip, the correspond- 
ing renormalized operator is given by 



(26) 



vergences by taking fi ^ v. 

ii) In presence of a charm- up GIM mechanism, the 
renormalized operators have the form 0,Ea] 



6±(/i) = z±(/i) 



+ (mc ~ mu)(ms - md) 5750! 

a 

nC|(5o) 



(27) 



where the Oi are operators of dimension six, the 
mixing of which is induced by the explicit chiral 
symmetry breaking of the lattice action and have 
been discussed in sec. |[ In this case, it has been 
impossible so far to get the matrix elements of the 
renormalized operators with resonable accuracy. 

iii) By defining 5*0 = hyDghy and K = hy{i D 
)'^hy, the renormalized operators have the form 

Soip) = Zhifl){hyD„hy + ^jMlhyhy^ 

k{p) = ZD2{p)(hy{iDfhy (28) 



In this case it is easy to get rid of the power di- 



The subtraction procedure necessary to obtain a 
finite 5*0, which is related to the definition of the 
renormalized quark mass, is straightforward [ p8[ . 
An accurate determination of the quark mass, 
however, requires the control of high orders in 
perturbation theory |39| ; the definition of a finite 
K, necessary to obtain the heavy quark kinetic 
energy is much harder and a small error on the 
final answer difficult to obtain | |3^ . 
The use of perturbation theory in the calcula- 
tion of the mixing coefficients necessary to sub- 
ract power divergencies gives unreliable results. 
A detailed discussion, with several examples can 
be found in |3|] . For the calculation of logarith- 
micaly divergent or finite mixing coefficients in- 
stead, perturbation theory can be used, at least 
in principle. In practice, however, due to large 
higher-order corrections present in the expansion 
of the lattice perturbative series, the results are 
not very accurate and it is for this reason that 
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the non-perturbative methods discussed in sec. ^ 
have been developed. 

All the problems connected to the construction of 
finite operators, out of the bare (power) divergent 
lattice ones, can be overcome, at least in princi- 
ple, by using the OPEwOPE. I start by recalling a 
result presented at the Latt97 , and then dis- 
cuss new proposals, and feasibility studies, which 
appeared this year. 

Weak Hamiltonian The first example of OPE- 
wOPE is related to non-leptonic weak decays. 
The standard construction of the weak Hamilto- 
nian begins with the expression 



4^ I d'xDZ{x;M„)T U{x)jU0) 



(29) 



where D^{x; M^) is the VF-boson propagator 
and JpL is the (left-handed) hadronic weak cur- 
rent. One then performs the OPE of the product 
of the two currents in eq. (p9|) 



{B\nZ\A) 

Gp 

71 



;b|6«(m)|a) (30) 



where di is the dimension of the operator O*^*-' (/i), 
and the functions Ci{fi, M^^,) result from the in- 
tegration of the Wilson expansion coefficients, 
Ci{x;fi) (defined in eq. ( ^2|) below), with the W- 
propagator. Schematically, suppressing Lorentz 
indices, one has 



(31) 



The sum in the expansion ( pO[ ) is over operators 
of increasing dimension. In the following only op- 
erators with dimensions di < 6 are considered, 
since the contribution from operators with c?i > 6 
is suppressed by powers of 

All the intricacies of operator mixing in the def- 
inition of the finite and renormalized operators, 
O*^'^ (/i), come about because the integral in ( p9| ) is 
extended down to the region of extremely small 
X. The complicated mixing for the 6(*)(Ai)'s m 
terms of bare operators arises from contact terms 
when the separation of the two currents goes to 
zero (i.e. when |x| is of the order of a). This 
observation suggests that a simple way to avoid 



these complications is to define the renormalized 
operators by enforcing the OPE for distances \x\ 
much larger than the lattice spacing a. This is 
done by computing directly the T-product of the 
two currents rather than the matrix elements of 
the local operators on the r.h.s. of eq. (|30|). For 
this reason I call this method OPEwOPE. We 
imagine proceeding in the following way: 1) Take 
the T-product of two properly normalized weak 
currents, JpL(x) J^^(O). If required these currents 
can be improved. 2) Measure the hadronic ma- 
trix element {B\T[JpL{x)jl^{0)\\A) in a Monte 
Carlo simulation, as a function of x for |a;| — > in 
the region a <C |a;| <C ^qcd- 3) Extract the num- 
bers (Slot*) {fi)\A) by fitting in the region defined 
above the x-behaviour of {B\T[JpL{x)jl^{Q)]\A) 
to the formula 



(BIT 



MJIlW \A) 



c,{x;fi){B\0^'^\^i)\A), 



(32) 



where the Wilson coefficients Ci{x; fi) are deter- 
mined by continuum perturbation theory using 
any standard renormalization scheme. Since we 
only consider operators of dimension 6 or lower, 
the T-product differs from the right-hand side of 
eq. ( |3^ ) by terms of OdxpAgp^,), which is an es- 
timate of the size of the systematic errors in this 
procedure. 4) Insert the numbers {B\0^'''> {n)\A) 
determined in this way into the expression for the 
matrix elements of Ti^ in eq. ( |30| ) . 
For the implementation of this procedure, what 
is required is the existence of a window, in which 
the distance between the two currents is small 
enough, so that perturbation theory can be used 
to determine the expected form of the OPE, but 
large enough that lattice artifacts are small. The 
existence of a window is a necessary conditions in 
all applications of the OPEwOPE. 
The method determines directly the "physical" 
matrix elements of the operators appearing in 
the OPE of the two currents, i.e. the matrix 
elements of the finite, renormalized operators 
0''^(/i), without any reference to the magnitude 
of the VF-mass. Thus we do not need to probe dis- 
tances of 0{1/M^^.) with lattice calculations. The 
/i-dependence of the matrix elements of the oper- 
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ators O^^^ (n) is given trivially by that of the (per- 
turbative) Wilson coefficients, Ci{x;^). It com- 
pensates the related /i-dependence of the func- 
tions Ci{y,,M^.) in such a way that the l.h.s of 
eq. ( ^0|) is independent of the choice of the sub- 
traction point. A similar comment holds for the 
dependence on renormalization scheme. 
Non-linear cr-model A feasibility study of this 
method, in a simplified case, has be presented 
by S. Caracciolo in a parallel session at this 
Conference They have considered the two- 

dimensional non-linear cr-model, for which it is 
possible to define the following conserved currents 

jf{x) = i {a%x)d,a\x) - a\x)d,,a\x)) (33) 

where a = (cr^, . . . , ) with the condition <j{x) ■ 
(j(x) — 1. The relevant expression is 

(p| Ej^'''(^)^'^"''(o)b) = CT{x-y^)(p\f^M (34) 

where \p) is some external state of momentum p 
and T^^l, the renormalized energy-momentum ten- 
sor (T^y — d^u • d^a). This operator has zero 
anomalous dimension, but it is subject to a finite 
renormalization (analogous to Zy a in the case of 
Wilson fermions) on the lattice. 
In the OPEwOPE approach, one computes non- 
perturbatively the matrix elements of the T- 
product of the original currents, corresponding 
to the l.h.s of (^) and the Wilson coefficients (in 
this case there is only one operator correspond- 
ing to Ct{x; n)) in perturbation theory. From 
the knowledge of these quantities, one can then 
extract the desired matrix element of the local 
operators, i.e. {p\T^y\p). In the example consid- 
ered here, we may also compute directly the ma- 
trix element of T^^ and check to what extent the 
procedure works. In [^0| they found a large win- 
dow of distances for which the OPE is verified, 
when continuum perturbation theory is used in 
the calculation of the Wilson coefficient, giving 
the correct matrix element independently of the 
external states. 

The quark propagator The simplest T-product 
that can be considered is the quark propagator in 
a fixed gauge, for example the Landau gauge. Let 



us consider for example the following quantity 

I]2(x) = (0|V'(x)S^(0)S|0), (35) 

where a and A are spin and color indices respec- 
tively. At short distances, i.e. as |x| — > 0, we 
have 



T;2{x) C^{x]ii){ip^{n)) +C„l{x■,^l) 



-o 



(36) 



where C^(a;;/i) and Cmix; are the Wilson 
coefficients, {r/j'ip{fj,)) the quark condensate and 
©(a/jxl^) represent (gauge non-invariant) lattice 
artifacts [^. The Wilson coefficients, which are 
gauge-dependent, can be computed in any con- 
tinuum renormalization scheme, for example in 
the MS scheme |42]. In the chiral limit, by 
fitting S2(x) to (36), and using C^{x]ijl) and 
C'm{x; fi) from perturbation theory, one can ex- 
tract the value of the condensate and the renor- 
malized quark mass in the same scheme as the 
Wilson coefficients. It would be very interesting 
to compare the values obtained in this way, with 
those obtained with other techniques. Note that 
if we use the Fourier transform of S2, there are 
further contributions from contact terms, which 
can make the extraction of the condensate and of 
the quark mass more complicated. 
The shape function The knowledge of the 
shape function /(fc+) is a fundamental ingredi- 
ent for the extraction of \Vub\ from the end-point 
of the lepton spectrum. The same function also 
enters the calculation of the photon spectrum in 
radiative B decays. To give an explicit exam- 
ple, the differential distribution in semileptonic 
decays is given by 



dV 
dE'i 



Mb 



dT 



PM 



dEf 



dk+ f(k+) 



(37) 



where ml — Mb 



k+. 



The physical rate can be derived from the imag- 
inary part of the forward matrix element of the 
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T-prodiict of two weak currents 
1_ 



-Imi 



B(v)\T \x) rifi) \B{v)) 



(38) 



In the standard approach, one apphes the OPE 
to the T-product above, and W^'^ is written in 
terms of the following matrix elements of local 
operators 

Mn (X {B{v)%-i%iD'^') . . . {iD'^-)K\B{v)) (39) 

which correspond to the moments of the shape 
function 

I- Mb 

Mn = / dk+klf{k+) (40) 

With the OPEwOPE, it is possible to obtain the 
full shape function. The final result of |43| is very 
simple. One study the ratio of two correlation 
functions 



(41) 



lim 



-Mb* 



where the numerator corresponds to insertion of 
the three-dimensional Fourier transform of the 
(Euclidean) T-product for t > 

and the denominator is the usual S-meson prop- 
agator 

stfM = mlBM'^PB=o{~um (42) 

^PB (0 is the B interpolating field with definite 
spatial momentum pB 



(43) 



^Psit) = J rf'a;e-^P--^$B(x,i) 
In terms of f{k^), we have 



(44) 



1 

2^ 



[k++^ 



dk+ f{k+)- 



where Q'^ = k^ 



kl + Q^ 



k+ + jQ^. By 



a suitable choice of the Lorentz components n 
and of the currents and of the spatial momen- 
tum Q, and by studying the time dependence of 
W^^lt, Q), we can unfold the integral above and 
extract the shape function. As explained in ]4^ , 
the same method can be used to extract the struc- 
ture functions of deep inelastic scattering. 
Light-cone wave-funtions The light-cone wave 
functions allow to predict form factors relevant in 
many exclusive processes, such as electromagnetic 
elastic scattering at large momentum transfer, or 
exclusive semi-leptonic decays as B ^ w {B ^ p) 
and B K*-f decays [|4)-@. The knowledge 
of the relevant light-cone functions would allow 
the determination of the form factors at ^ 0, 
a kinematical region which is not accessible with 
standard lattice techniques. In |^^, it was shown 
that they can be computed, analogously to the 
shape function, by suitable ratios of three- and 
two-point correlation functions of Euclidean T- 
product. In this case the three-point function cor- 
responds to the insertion of the three-dimensional 
Fourier transform of the T-product 

(^-|r[(i»exp*.''o 7^75^(0)1 10) (45) 



where = (g^ +p^)/ y— and 



The 



OPE of the T-product can be written in terms of 
matrix elements of local operators 

Mnoc (7r-|J(0)7^75(iT'^^)...(iT'^")w(0)|0) (46) 

which, as for the shape function considered be- 
fore, corresponds to the moments of the ligh-cone 
wave-fuction 



du m"$^(-u) 



(47) 



The matrix elements of the local operators cor- 
responding to the first two non-trivial moments 
were studied, using the standard approach in . 
With the OPEwOPE one can obtain the fuU fight- 
cone wave function $ without the problems asso- 
ciated with the renormalization of the local oper- 
ators of eq. (|4€ 
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